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Abstract 

We apply the principles discussed in an earlier paper to the construction 
of discrete time field theories. We derive the discrete time field equations of 
motion and Noether's theorem and apply them to the Schrddinger equation to 
illustrate the methodology. Stationary solutions to the discrete time Schrddinger 
wave equation are found to be identical to standard energy eigenvalue solutions 
except for a fundamental limit on the energy. Then we apply the formalism 
to the free neutral Klein Gordon system, deriving the equations of motion and 
conserved quantities such as the linear momentum and angular momentum. 
We show that there is an upper bound on the magnitude of linear momentum 
for physical particle-like solutions. We extend the formalism to the charged 
scalar field coupled to Maxwell's electrodynamics in a gauge invariant way. We 
apply the formalism to include the Maxwell and Dirac fields, setting the scene 
for second quantisation of discrete time mechanics and discrete time Quantum 
Electrodynamics. 

1 Introduction 

THIS paper follows an earlier paper on the principles of discrete time mechanics as 
formulated for point particle theories, hereafter referred to as Paper I |L|. These 
principles are applied in this paper to systems described by fields (p a (t, x) in one time 
and 3 spatial dimensions, where a denotes spin component and field types. Units will 
normally be taken to give c = h = 1. The summation convention will apply only to 
repeated small Greek indices except where otherwise stated. 

In discrete time mechanics dynamical variables take on values at times t n , n = 
0, 1, 2, ...N, with the discrete time intervals At n = t n+ i — t n , n = 0, 1, 2, ... given by 
At n = T , where T > is the fundamental time unit of the system. An important 
principle of our formulation is that the mechanics is regarded as exact regardless of 
the magnitude of T, and not an exercise in approximation. This means that invariants 
of the motion are constructed to be conserved precisely and not up to some given 
order in T. 



As in Paper I our discrete time approach takes its cues from continuous time 
mechanics, but ultimately is more than an approximation to such a formalism. In 
general the continuous time Lagrangians L for the systems of interest to us in this 
paper are the spatial integrals of Lagrange densities C of the form C = C((f, d t <£, Vy>). 

Important concepts in our methodology are those of the system function and 
virtual path. In our formalism the system function has the status occupied by the 
Lagrangian in continuous time mechanics. We may calculate equations of motion and 
conserved quantities if we know the system function. We may also use the system 
function to define conjugate momenta and to quantise our theory. We construct 
the system function from the continuous time Lagrangian by integrating over certain 
trajectories, which we refer to as virtual paths. These are not solutions to the normal 
Euler-Lagrange equations of motion but paths chosen according to specific criteria 
such as gauge invar iance. 

In our approach we shall discuss discrete time mechanics based on a temporal 
lattice with successive points given by 

tn+l —t n — T, (1) 

where T is non-zero and independent of n. If tp (t) is a dynamical variable depending 
on time we will normally define a virtual path ip from ip n = <p (t n ) to p n+ i = (p (t n +i) 
in <y?-value space by the linear rule 

<f (t) = Vn+l + V n , (2) 

where A runs from to 1 and A = 1 — A. This means that over the interval [t n , £ n +i] 
the time t is related to A by the rule 

t = t n + XT = \t n+ i + \t n . (3) 

This prescription is modified when we include gauge invariance, as discussed in §8. 
With the above virtual path prescription time derivatives are naturally turned into 
differences, i.e. 

^ (i) = ^±LZ^. (4) 

It is possible to use a different temporal lattice. A popular choice in recent years 
has been the so-called g-lattice, where we define t n+ i = qt n , where q is real and not 
equal to unity. Then we would define ip n = p> (t n ) = p> (q n t ) , t ^ 0. The virtual 
path prescription on the interval [t n ,t n+ i] would be the same as (Q) with 

t = t n + X (t n+ i — t n ) = Xt n+ i + At n , (5) 

and so in the interval [t n , t n +i] the temporal derivative would be replaced by the 
so-called g-derivative 

(f) = ¥n+l ~ = V (fin) ~ <fi (t n ) 

dV { ' t n+1 -t n (q-l)t n 1 ' 
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We may go further and accommodate the so-called symmetric derivative given by 

9 ~ ,.x _ y? (qt n ) - if (q-H n ) , . 

ai v{t)= («-«-')*. ' (7) 

but in this particular case the virtual temporal path (viewed in terms of successive 
intervals) is not continuous everywhere. We may even generalise in the manner of 
Klimek || and define t n +\ = <p (t n ), where <ft is some chosen real valued function and 
then the above prescription gives 

dP {t) - Htn)-t n ■ 

The regular temporal lattice is given by the function <p(t) = t + T in Klimek's 
formulation. We shall use this lattice throughout this paper, applying it to the 
Schrodinger wave equation, the free neutral scalar field, the free charged scalar field, 
the Maxwell fields and potentials and to the Dirac field, and show how to couple 
these fields in a gauge invariant way. The quantisation of our discrete time field 
theories will be discussed in the next paper of this series, Paper III, following the 
principles discussed in Paper I. A detailed discussion of quantum electrodynamics 
will be presented in Paper IV of this series. 

The plan of this paper is as follows. First we review Hamilton's principal function 
in continuous time field theory and then discuss the related concept of system function 
in discrete time field theory. We derive the field Cadzow's equations of motion and 
discuss the Maeda-Noether theorem for the construction of invariants of the motion. 
Then we apply our formalism to the Schrodinger equation, as an example of a T>Q 
process, where Q represents the process of quantisation and V represents the process 
of discretisation of time. The resulting theory is a classical field theory which is not 
the same as the result of the QD process discussed in Paper I. 

Then we turn to relativistic theories. The T> process breaks Lorentz invariance, 
but in a less emphatic fashion than space-time lattice theories such as those discussed 
by Yamamoto et al. || || . The systems we study here are the free neutral and charged 
Klein-Gordon fields, the Maxwell potentials and their coupling to the charged scalar 
field, and the Dirac equation. This sets the scene for second quantisation, i.e. the 
QVQ process, which is discussed in Paper III. 



2 Hamilton's principal function 

In this section we review some aspects of continuous time mechanics field theory 
which are relevant to our method of discretisation. Following the principles outlined in 
Paper I we shall consider a dynamical field system for which we have knowledge of the 
dynamical variables <p a at a sequence of times t n , n = 0, 1, 2, ...N, where t n+ i = t n + 
T. An important construction in continuous time mechanics is Hamilton's principal 
function S n (T) , defined as the integral 

tn+l 

S n (T) = J dtfcPxC (9) 
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over the classical solution to the field equations which satisfies the boundary condi- 
tions 

< (x) = cp a (t n , x) , < +1 (x) = cp a (t n+1 , x) . (10) 

S n (T) is a functional of <^> n and <p n+ i and is the continuous time analogue of the 

system function in discrete time mechanics. If the momentum n a conjugate to ip a is 
ac 

d(p a 



defined by n a = and the Hamiltonian H defined by 



H = J d 3 x{n Q ^ a - £} (11) 
then for an infinitesimal variation 5ip a of the fields we find 



SS n (T) 



-5tH + J d 3 xn a 5(p a 



tn+l 



(12) 



where the symbol = denotes an equality holding over the true or dynamical trajectory. 
From this we deduce the functional derivatives 

s. on ton 
^n ( x ) ^n+l ( x ) 

and the partial derivatives 

Q an f) on 

^_ = #n = (14) 

Ot n Ot n+ i 

where H n is the value of the Hamiltonian at time t n . 

By the process of constructing S n for each of the intervals [t n t n +i] , n = 0, 1, N— 
1, we arrive at an action integral from t to tjy which is equivalent to the action sum 

N-l 

A N =J2 s n - ( 15 ) 

n=0 

If the end point field values are held fixed (at times to, t^) then this sum depends on 
the field values tp^ (x) at the intermediate time t n , n = 1,2, N — 1. If now we vary 
these intermediate fields and apply a variational principle we find the equations 

IS™" 1 + S n \ = 0, < n < N. (16) 

6<f% (x) I / c 

These are the continuous time mechanics analogues of Cadzow's equations of motion 
f| in discrete time field theory, discussed in the next section. The above equations 
([16|) essentially ensure the continuity of the conjugate momenta at the end points of 
the sub-intervals. 
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3 Discrete time field theory 



We turn now to the discretisation of Lagrange based field theories based on the 
construction of the system function F n using the principles outlined in Paper I. By 
definition, the system function is the integral 

rtn+i r „ 

F n = dt d 3 x£(<p,d t ¥,Vip) (17) 

Jtn J 

with some choice of virtual path (p in field value space. Once the temporal integral has 
been carried out the system function is a function of the field values v?n( x ) = <P a (tn, x ) 
at times t n , t n+ i and their spatial derivatives at those times. 
Given the system function, the action sum A N is given by 

N-l 

a n = J- F n (18) 

n=0 

and it is this which replaces the action sum (|T3J) in continuous time mechanics. We 
derive the equations of motion using flTB] ) and a variational principle in the stan- 
dard fashion || . Consider infinitesimal variations of the fields with fixed end-points, 
5<£ (x) =5cp N ('x) = 0. Assuming the fields fall off rapidly at spatial infinity, we apply 
Cadzow's action principle to obtain the discrete time field equations of motion 

{F n + F n - 1 } = 0, 0<n<N, (19) 



6(ft (x 
which are equivalent to 

where T n = J-((p n , Vy n , (fn+u V^n+i) is the system function density defined by 

rtn+i 

T n = \ dt£(g>,d t &vg>). (21) 

The construction of constants of the motion is straightforward using Noether's the- 
orem applied to discrete systems. Consider an infinitesimal transformation 5cp n , S(p n+1 
of the fields which leaves the system function unchanged. Then using the equation 
of motion we deduce that the quantity 

C n = /d 3 x{^-V-U^)(5< (22) 

is conserved on dynamical trajectories. In Paper I we referred to this as the Maeda- 
Noether theorem |J. 

For linear momentum consider a transformation of the fields by the infinitesimal 
shift 

<(x) -^(x) = <(x) + 5a-V<(x). (23) 
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Assuming that the Lagrange density is not explicitly dependent on position then we 
find that the conserved linear momentum P™ is given by 




For orbital angular momentum, consider a transformation of the fields by an 
infinitesimal rotation 

y^(x) -^(x) = <(x) + (Slj x x) -Vy^(x) + 5u> ■ 4#<pP n (x) , (25) 

where s a/3 is the transformation matrix for the spin indices. From this we may con- 
struct the conserved angular momentum 

L»,/ d 3 x(xxV< + S ^){^-V.(^)}. (26) 

Other conserved quantities such as electric charge are just as readily obtained by 
the same method. 



4 The discrete time Schrodinger equation 

In Paper I we discussed the possible non-commutativity of the processes of quantisa- 
tion Q and discretisation V. In that paper we presented an approach to Q after the 
application of D to point particle Lagrangians. Here we shall consider applying our T> 
to a case where Q has already been carried out. The first application of our approach 
to the temporal discretisation of field theories is chosen to be the Schrodinger wave 
equation for a particle in an external time-independent potential. 

First we note that the Schrodinger equation may be derived as a classical field 
equation in continuous time mechanics using the wave function \I/ and its complex 
conjugate as dynamical field variables with the Lagrangian density 

*2 ^_ _^ 

C = ph(V*dtV - d t ^*y) * V • V * - V (x) . (27) 

The Euler-Lagrange equation of motion 

PC 8C dC 

leads to the usual Schrodinger equation 

ihdtV = + V (x) * (29) 

c 2m 

and similarly for the complex conjugate wave-function. 

Turning to our discretisation process, we define our virtual paths by 

tf n (f,x) = A^ n+1 (x) + (x) , 

*;(f,x) = a*; +1 (x) + a*;(x), 

f(A) = t n + \ (t n+1 - t n ) , (30) 
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for t G [t n ,t n+ i], where \l/ n (x) = \& (i n ,x) , etc. and < A < 1. Then the system 
function turns out to be 

2 I n n+1 ^n+l^nj 

-| {2 |^ n+1 | 2 + 2 |^ n | 2 + % +1 v n + K^n+i} v (x) 

{2 |W n+ i| 2 + 2 |W n | 2 + K+i V ■ V^ n + *;V ■ V^ n+ i} .(31) 
Using (pT]), Cadzow's equation of motion for \1/ is 

«<*"+■ -*-'> = f-^ + V (x)l (W4tt„ + *„-.) (32) 
2T c [ 2m j 6 

and similarly for the complex conjugate. This is our discretisation of the Schrodinger 
equation and represents the process T>Q applied to a standard continuous time me- 
chanics point particle system. 

Although in this theory there is no direct concept of a Hamiltonian, this does not 
mean that there are no conserved quantities analogous to continuous time energy. 
We define a stationary state solution to the discrete time Schrodinger equation (|32|) 
to be of the form 

n (x) = e- itnT/n iP (x) (33) 
where e may be thought of as a discrete time energy. Then we find (|32|) reduces to 

^(x) = |-^ + V(x)|^(x) (34) 

where the eigenvalue E is given by 

3h sin (eT/h) 
T{cos(eT/ft) + 2}' 1 ' 

A plot of y = ET/h as a function of x = eT/h shows a periodic graph which passes 
through e = with a slope of unity, so that for low energies (compared to h/T) e is 
virtually identical with E. Then the graph rises to a local maximum when 

e = — (36) 
3T v ; 

at which point E has the value E max = :y ^. If, as we imagine, T is a very small 
timescale, such as the Planck time, then -E max is in practical terms far beyond any 
energy scale encountered in the laboratory or in any astrophysical or cosmological 
context, except possibly close to the Big Bang. Beyond the local maximum at e = ^ 
the graph falls to zero and goes negative, followed by a local minimum and a return 
to zero. Thereafter it repeats this overall pattern. We note here the possibility of 
resonances, where the values of e differ by multiples of These vanish in the limit 
T -> 0. 
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Assuming < E < and disregarding the resonances, there are two distinct 
values of e which would give the same value of E. Of these, the lower value would 
correspond to a normal physical state energy in practice whereas the higher value 
would be associated with a state with an enormous temporal oscillation factor. We 
expect that such a state would not be created in the laboratory under normal situa- 
tions, because the scattering amplitudes to go to such states oscillate very rapidly in 
time and this would lead to cancellations to zero in detailed scattering calculations. 

As far as the resonances go, these give stationary solutions which differ by at most 
a sign to one of the two mentioned above, and so do not represent new states and 
can be ignored. 

An interesting question is whether every eigenfunction of the continuous time 
equation (|4]) generates a stationary state solution to (^) of the form (|33|). The 
answer is no if we require e to be real. For energy eigenvalues greater than 
solutions to (|35D will develop an imaginary component, giving discrete time wave- 
functions which grow or decay exponentially, and such wave-functions cannot be 
regarded as being stationary. An analogous cut-off phenomenon will be seen when 
we discuss the Klein-Gordon equation in §5. 

To find a conserved charge we first rewrite the equation of motion ([32]) in the form 



{3ih - T5)f n+1 = (3ih + TsQ ^ n _! + 4T sV n 
1 (-3ih - T*S) = (-3ih + T*S ) + 4TV* n *S (37) 



where S = ^-V 2 + V (x) . Now consider the global gauge transformation 

#„->e*tt n , %^e~ i6 ^ n (38) 

and apply the Maeda-Noether theorem discussed in §3. Then we find the invariant 
of the motion 

QT = j rf 3 x + l 0)^ n+l + K + i (~ - ^g) , (39) 

which reduces to the standard total probability in the limit T — ► 0. With the equations 
of motion ( j37D we find 

or = Q n - 1 . (40) 



From the charge (|39D we construct the charge and current densities, given by 

Pn = *;(i + f^)*n + l+K +1 (}- *Sjx)*n 

k = ^{^V^ + ^V^ + KV^-!}, (41) 

which satisfy the discrete time equation of continuity 

Pn ~ T Pn - 1 + V j n = 0. (42) 
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The density p n above is not positive definite except in the limit T — > 0, when the 
normal density is recovered. This is analogous to what happens in continuous 
time mechanics to the charge density for the charged Klein-Gordon equation before 
we take the limit c — > oo. One of the reasons occasionally cited for Schrodinger's 
rejection of a relativistic wave equation in 1925-26 is because of just this point, there 
being no evidence at that time for positive electron charge densities. We imagine 
that if discrete time mechanics had been the accepted classical mechanical paradigm 
prior to wave mechanics then the Born probability interpretation of Schrodinger wave 
mechanics would probably not have been proposed. 

If we wish to couple electromagnetic fields in a gauge invariant way to the Schrodinger 
equation we need to choose a modified virtual path. Such a path is used for the 
charged Klein-Gordon equation discussed below, and so this is left as an exercise in 
the case of the Schrodinger equation. 

5 The Klein-Gordon field 

5.1 The continuous time system 

An important system is the free scalar field system with Lagrange density 

C = \d^d»ip - i/xV, (43) 

where we shall take c — h — 1 and three spatial dimensions. Hamilton's principal 
function is given by 

S n = 1 f d 3 p Po I 3 3 ip . (x _ y) (44) 

c 2 J ( 2vr ) 3 sin (p T) J y y 1 

{ [fn+1 ( X ) ¥n+l (y) + Vn ( X ) Vn (y)] C °S (j> T) ~ 2y?„ +1 (x) Vn (y) } , 

where po = VP ' P+A* 2 - Now compare this principal function with that for the con- 
tinuous time mechanics harmonic oscillator, where the phenomenon of recurrence 
occurs, as discussed in Paper I. We note that in the field theory, recurrence occurs 
in momentum space. By inspection of (|44|) , singularities in the momentum space 
integrand appear to arise whenever poT = rmc, m = 1, 2, ... To understand this here 
we first construct the conjugate momenta in terms of the end-point field values using 
(O) . From the principal function ([01) we find 



7T, 



1 = / ■ 7 ° / rf V e ip - (x - y) [p n+1 (y) - (y) cos (p T)} , (45) 

J (27r) sin (poT) J 1 J 

which satisfies the expected relation 

lim7r n (x) = v? n (x) (46) 

if we assume the fields vary smoothly with time. 
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Important constructions are the momentum fields 

A(p) = J d 3 xe ipx {p oV + nr} (47) 

A*(p) = Jd^e-^ipoip-in}, (48) 

which are time independent modulo the Klein-Gordon equation. Using we find 

An (p) = — / rf 3 x e ^oT- ip .x (x) _ (x) } (4Q) 

and similarly for its complex conjugate. The time independence of the above mo- 
mentum fields means that 

A n (p) = A n+1 (p) . (50) 



Moreover, the momentum fields are not singular, as is evident from (47, 4q) , so that 
in the limit poT — > mir the apparent singularity in A n (p) must be cancelled by 
recurrence behaviour in the Fourier transforms 

ip n (p) = f d 3 xe-*->„ (x) . (51) 



We deduce from ( p^Oj ) that the Fourier transformed fields y? n (p) satisfy the equation 
of motion 

Vn+i (P) + <P n -i (P) = 2 cos (PoT) y n (p) . (52) 

This shows that the continuous time Klein-Gordon field behaves like the discrete 
time harmonic oscillator discussed in Paper I and the discrete time Klein- Gordon 
field discussed next, but with one important difference. The magnitude of the factor 
cos (poT) in the above equation of motion never exceeds unity and this ensures that 
the motion is never hyperbolic. There is no need therefore for a momentum cut-off 
in the continuous time theory. 

5.2 The discrete time Klein-Gordon field 

We now consider discretising the continuous time system with Lagrange density ( f43|) . 
The virtual paths are of the form 

^(x) = A^ n+1 (x)+A^(x), 0<A<1, (53) 

where A = 1 — A. With this choice of virtual path the time derivative d t may be 
replaced by the operator T~ 1 d\, which acts on up as a difference operator, viz 

^ = ^±lf^. (54) 
The system function F n = F (tp n , <p n+1 , Vy2 n , Vv? n+ i) is then the spatial integral 

F n = J d 3 xr(x) (55) 
10 



of the system function density jF n (x), given by 

i 

J- n (x) =T J d\C(<p,dp$). (56) 



We find 



{iPn+ \ T ^ - f (|v^ n+1 | 2 + V Vn+l ■ V Vn + |V^| 2 ) 

/i 2 T 



(<Pn+l + Vn+l^n + ) (57) 
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from which Cadzow's equation of motion (|20|) gives 



^2 ^ g c V / 

where K = —V 2 + /i 2 . This equation can be written in the form 

D<p n+1 + AD Vn + ~Dip n _ 1 = -ip n (59) 

where 

-> 6 + T 2 K , . 

s -in^' (60) 

which is useful for proving various quantities are conserved. For example, the linear 
momentum obtained using ([24]) is found to be 

j d^ n VD Vn+l (61) 

and using (|59D we find 

pn = pn-1 ^ 



c 



Likewise, the orbital angular momentum is conserved and given by 

L" = - f rf 3 x (x x Vy? n ) ^ n+1 . (63) 

These invariants of the motion exist because we have not destroyed Euclidean in- 
variance, that is, there is still spatial translational and rotational invariance in our 
approach. This will occur in all our discrete time field models based on special rela- 
tivistic Lagrangians. We note that a similar equation was found for the scalar field 
in the k— Poincare theory of Lukierski et al. ||, with the important differences that 
their lattice parameter corresponding to our T is imaginary and that they are dealing 
with continuous time throughout. 

We turn now to particle-like solutions to the equation of motion. Consider the 
Fourier transform 

£ B (p) = /^e" ,p \(x). (64) 
11 



The equation of motion (|59|) becomes 
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/3£> n+ i + 4/3£ n + /3^ n _ x = -£> n , (65) 

where (3 = 6+ gy £2 and i£ = \/// 2 + p • p. Now define the momentum functions 

a n (p) = if3e me [<p n+1 - e i0 cp n ] , (66) 
<(p) = _^ e - me - e"^;] , (67) 

where cos# = i], sin 6 = y/1 — rj 2 with t] given by the ratio 

a 6-2T 2 £ 2 

V — Oi = — . (68) 

' (3 6T K ' 



These are closely related to the continuous time momentum functions ([49[ ) and to 
the ladder operators discussed in Paper I when we identify e %e = \i = r] + — rj 1 . 
Particle states in discrete time field theory correspond to elliptic type wave behaviour, 
which requires < 1. Otherwise, physically unacceptable hyperbolic behaviour 
occurs, as discussed in Paper /. This is equivalent to the condition 

TE < 2^3, (69) 

i.e. 

P"P + /i 2 c 2 <^ (70) 

in Standard International units. The fundamental conclusion from this analysis is 
that energy and momentum are bounded above for physical particle states in our 
approach. In other words, in discrete time field theory, the spectrum of acceptable 
particle states has a natural cut-off which can be made as large enough to avoid a 
violation of observed particle data as necessary by choosing the time interval T small 
enough. 

An interesting point is that for a given T, there can be no scalar particle species 
with rest mass greater than 2\^3h/Tc 2 , but this cannot be regarded as of interest to 
experimentalists at this time. 

We note that by construction the momentum fields are time independent and 
satisfy 

a n (p) = On-i (p) , a* (p) = <_! (p) . (71) 
This allows us to construct a Logan invariant of the motion 

_fp_ (12 

(2tt) 3 It 2 

which in the limit T — > becomes the standard field theory Hamiltonian, using the 
result 

lira a n {p)=i f d 3 xe iEt ~ ip ' x [cp(t, -x)—iE(p (t, x)] = A(p) . (73) 

1 — >U, nl — >t J 

for those trajectories for which the limit exists. It is important to understand however 
that there is no notion of a Hamiltonian in our approach, since there is no possibility 
of infinitesimal translation in time, except in the above limit. 



C n = ™-A* -P-P K(PK(P) (72) 
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6 The free charged Klein Gordon equation 



For this system we consider the continuous time Lagrange density 

C = d^d^if — fi 2 (p*(p. 
The virtual paths are given as before by 

<p = Vn+i + v„, v* = Vn+i + v;, < A < 1, 



(74) 



(75) 



where we suppress the dependence on x. Then the system function density is found 
to be 



T 



// 2 T 
6 



T -12 

6 



+ 2 |^| + Vn+lVn + VnVn+l f ■ 



(76) 



From this Cadzow's equations of motion give 
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which can be written in the form 




0. 



£>^ n+1 + 4L>y? n + Dip n _ x 



T 



D Vn+l + 4D(f* n + D {p* n _ 1 = —iPn, 



(77) 
(78) 

(79) 
(80) 



where D is given by fl6"0"|) . 

The linear and angular momenta are easy to construct so we turn to the new 
feature, global gauge invariance. The system function density ([76]) is invariant to a 
global gauge transformation of the fields, i.e. 

<p n -Kfif n = e~ l ^ n , ifl - ipt = e>;, (81) 

where x is independent of time and space, so using the Maeda-Noether theorem 
discussed in §3 we find the conserved charge 

Q n = * j 'd 3 x{v* n D Vn+1 -v* n+1 D Vn } (82) 

= ijd?x{ip* n *Dipn + i-<Pl + iDiPn}- (83) 
This is real and global gauge invariant. Using the equations of motion we find 

Q n = Q n -\ (84) 
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By inspection there are two possible candidates for a charge density, denoted by 
p*~) and p( + \ given by 



Pi ) = il P*n Dl Pn+l- il Pn+l Dl Pn, 
P ( n ] = UPn D( Pn+l- i( Pn+l D( Pn- 

These are related by a total divergence, that is 

iT. 



{+) = _(-) , ^ v 

l J n r n * (3 



(85) 
(86) 



(87) 



so that they give the same total charge Q n . Using the equations of motion we find 

:V • vV n + y?; V^ n „x + 4^; W 
-v- <Vy9„ +1 + y9; +1 vV„ + 4^; vV 



Hn rn—1 



p {+) - p (+ \ 



T 



^ 6 
^ 6 



(89) 



which are discrete time versions of the charge continuity equation, with corresponding 
charge currents given by 



Jn 



Jn 



—I 



(90) 
(91) 



7 Maxwell's equations 
7.1 Charge free equations 

We retain our unit system such that c = h = 1. Our discrete time formulation of 
charge free Maxwell's equations starts with the electromagnetic potentials and A 
which are used to construct the physical electric and magnetic fields E and B. A clear 
distinction has to be made here between the nature of the electric scalar potential 
(j) and the magnetic vector potential A. The former is associated with the temporal 
interval or link connecting times t n and t n+ i, whereas the latter is associated with 
the times themselves. This distinction also manifests itself in the difference between 
the physical electric field E and the magnetic field B, which are likewise associated 
with temporal links and endpoints respectively. 
Our definitions are as follows: 

The electric scalar potential associated with the link connecting time t n and t n+ i 
at spatial position x will be denoted by the symbol <p n (x) , rather than by (say) 
n+ i (x) . Although our notation suggests a bias towards t n at the expense of t n+ i, 

this is not really the case. The virtual path 4> n for the electric potential is defined by 



6 n (x) = n (x), 0<A<1. 



(92) 
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On the other hand the magnetic vector potential associated at time t n at spatial 
position x will be denoted by A n (x) and its associated virtual path A n is defined by 

A n (x) = AA n+ i (x) + AA n (x) , < A < 1. (93) 

A local gauge transformation involves the gauge functions \i which are associated 
with endpoints rather than links. A gauge function value at time t n and position x 
will be denoted by \ n ( x ) an d is assumed differentiable with respect to x. We shall 
see in the next section on charged scalar field electrodynamics that the virtual path 
for the gauge fields is given by 

X„(x) = Ax n+1 (x) + A Xn (x), 0<A<1. (94) 

A local gauge transformation is defined here by the replacements 

which reduce to the transformations 

0n - ^ = 0n+ W T ~ X " , (96) 

A n -> A; = A„,-Vx„. (97) 

Turning now to the physical fields, we define the gauge invariant electric and 
magnetic fields via the potentials: 

E n = -V0 n - ( A " +1 r ~ An) , B„eVx A n , (98) 
which give the homogeneous Maxwell equations 

V-B n = 0, V x E w + (B " +1 T ~ B " } = 0. (99) 

The gauge invariant system function density for the charge free system is defined 
as the integral 

i 

T n = T JdX (-ip n ^Fr) , (100) 



where F^ v = d^A^ — d u A£ are the components of the Faraday tensor. This gives 

= \K - I K +1 + B n+1 • B n + Bl) . (101) 
Applying Cadzow's equation to (|101|) we find 

V-E n = 0, (E "" r E "-' ) =Vx (B ° +1 + ^ + B "- l) (102) 
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Using ( p9|) and (|102| ) we find the physical electromagnetic fields satisfy the discrete 



time massless Klein-Gordon equations 

E w +i — 2E n + E n _i ^2 (E ra+ i + 4E ra + ~Ei n -i) _ n ,.„„, 

B n+ i — 2B W + B w _x 2 (B ra+1 + 4B„, + B n _x) _ , > 

7^2 v g 7 ^ ' 

Turning to the gauge potentials, we define the discrete time Lorentz gauge by the 
condition 

4> n ~ n _i , _ (An+i + 4A n + A„_i N 



T 6 
and then the gauge potentials also satisfy the discrete time massless Klein-Gordon 
equations 

A n+ i - 2A n + A n _i _ 2 (A w+ i + 4A W + A w _i) = 

^2 g c ' v > 



n+1 - 2 n + n _i ^2 



(0n+l + + 0n-l) 



/2 6 " ^ =0. (107) 

The total linear momentum for the free electromagnetic fields is found to be 

P" = J d 3 x{E„xB„ + ^;V5i +1 } . (108) 

Then using the equations of motion we find 

P" = P"" 1 (109) 
as expected. In the limit T — > the above expression reduces to the Poynting vector. 

7.2 Maxwell's equations in the presence of charges 

In the presence of electric charges the system function density ( |101|) is replaced by 

T n = \K-\ + Bn+i ■ B n + B^) - T<P nPn + TA n ■ j n (110) 

where p n (x) and j n (x) are the discrete time charge density and charge current re- 
spectively. The homogeneous equations remain unaltered but now the equations 
of motion become 

V-E n = Pn , Vx (Bn +1 +4B B + B w . 1 ) _ (E w -E w _ 1 ) =jwj (m) 

C (3 1 c 

These equations are consistent with the equation of continuity for electric charge 
given by 

Pn ~ T Pn ' 1 +V-j ra = 0. (112) 
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Finally, in a discrete time Lorentz gauge (|105| ) the potentials satisfy the equations 



A w+1 - 2A n + A n _x _ 2 (A w+ i + 4A ra + A w _i) _ . 

T 2 6 ~ Jn ' 

(j) n+1 - 2<f) n + n _! _ 2 (<t>n+l + 4< ^n + 4>n-l) _ (pn+1 + 4 Pn + Pn-1 



T 2 6 c 6 

(113) 

8 Scalar field electrodynamics 

We are now in a position to discuss the coupling of the Maxwell potentials to the 
charged Klein-Gordon field. In this case the virtual paths for the charged scalar 
field and its complex conjugate have to be modified so as to ensure that the system 
function is gauge invariant. We define 

~* = Ae -^„TA^* +i + Xe^ TX (p* n , A = l-A: 0<A<1. (114) 

where <ft n is the scalar potential discussed in the previous section. Then the gauge 
transformations are given by (|96|,|97D and by 



ipf n = e-*x*tp n , p>' n+ i = e- i,x " +1 ^i, (H5) 
which leads to the compact form 

(p f = e- iq *tp, ip*' = e iq *(p*. (116) 
Now define the gauge covariant derivatives 

D^=^ + iqA^ (117) 
and its complex conjugate. Under a gauge transformation we find 

ip^)' = e~ l ^D^. (118) 
Next we construct the gauge invariant Lagrange density 

= (D Q (p)* D Q (p- (D£)*-D£-//V£, (119) 

where 

/I \ P i( l<t> n T\ _ -iq<p n T\ 

Dot = (±d x + ^ q <P n ) (p = ^±i^ ^ (120) 



and 



tiip = (V - iqA)p 

= iqT\XS7(f) n {e* r V n+1 - e~^ TX Vn ) 

+ Xe iqKTX V(pn+i + \ e -i^nTX V(pn _ (121) 
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The system function is constructed as before from the integral 



F n = T J d\[ (p (p) * D V ~ ipv) * -D<^ - /i 2 (p*ip} 



(122) 



which gives 



T 
T 

6 

iqT 2 



2 T 3 



30 



-V0 n -V0 n U n <p 

n+l V^n 



n+l 



¥>n+l 



2 1 v^| 2 + ^:v^; +1 ■ v Vn + v< ■ v^ n+1 c/ r , 



+ 2 |^ n | 2 + ^">n+l^n + <Pn<Pn+l U r. 



U*<p* 
nrn 



12 \ -V0 n • (C/*V< +1 + V<) (i/ bVWi - ^ 

{ - <p* n ) V<p n ■ [(3A n+1 + 2A n ) £/ n ^ n+1 + (2A n+1 + 3A n ) cp r 

.)} 



+ - 

60 

+ [(3A n+1 + 2A n ) U*<p* n+1 + (2A n+1 + 3A n ) <p* n ] ■ V^ n (u n <p n+1 - <p, 
_tgT f (3A n+1 + A n ) • (< +1 V^ n+1 - V<^ 
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+ (3A n + A n+1 ) • «Vy?„ - V<p* n <p r 



+ (A n+1 + A n ) • [p* n Vip n+l U n - V(p* n (p n+1 U n + U*<p* n+1 V<p n - U*V<p* n+1 (p n ] } 
q 2 T \ (l2A2 +1 + 2A2 + 6A n -A n+1 ) ip n+l 



60 



+ (l2A2+2A2 +1 + 6A n -A n+1 )|^ n 



3A£ +1 + 3A* + 4A n • An+i) + <^>„,+i^n] } 



(123) 



where [7 n = e ig ^™ T '. This system function is gauge invariant. 

The application of the resulting Cadzow's equations to specific problems is left 
for possible future investigation. We make the following comments in passing. First, 
this system function looks complicated. However, that is a matter of notation, as the 
expression ( 122 ) is equivalent and much more compact. There is undoubtedly much 
dynamical content in this system function and calculations would no doubt require 
approximations to be made. However we recall here Gell-Mann's dictum that it is 
better to make approximations to exact equations than to solve exactly approximate 
equations. Also, most applications would themselves be better considered from a 
second quantised approach, which we shall study in Paper IV of this series. 

We may simplify the system function considerably by considering the special case 
when the external magnetic potential A n vanishes. Then the system function reduces 
to 



(I 



T 

T , 
— 2 



2^3 



30 
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+ 2 \ ( fn\ +U* ( P*n+l ( Pn + ( P*n ( Pn+l U n 



iqT 2 f (c/>; + i-<)V0 n -([/„V^ +1 + V^ n 

12 \ -V0„ • (u* v^; +1 + v<) (c/ n ^ n+1 - Vn ) ' 

which leads to the equation of motion 

U n(Pn+1 -2 (Pn + U*_ 1 i Pn _ 1 , /i 2 (C/^n + l+^n + ^-l^n-l) 



(124) 



T 2 



+ 



-V- 

g 2 T 2 

30 
iqT 



|(V0 n ) 2 [[/„^ n+ i - y? n ] - (V0 n _x) 2 [</?„ - C^_iy„-i] } 
{V0„ • (C/„V^ n+1 + V^J - V^.! • (V^ B + EdV^)} 



12 

V • { V0 n (tf n¥Wl - y?„) + V0 n _ x (y?„ - ICi^i) } • (125) 
In the limit T — > we recover the usual equation 

+ *g0) (9 t + *g0) y? + (V - V 2 ) y? = 0. (126) 

9 The Dirac equation 

We turn now to the continuous time Dirac equation 

iip = (-ia ■ V + /3rnJ ip, (127) 
which is obtained from the Lagrange density 

iip + cx ■ V ip — mip + (3ip. 



£ = -ijy itp + iot ■ V ip 



-iip + 



The virtual paths are taken to be 

and then we find the system function density 



(128) 



(129) 



T 
6 



{^l^n+l + ^+1 ^Vn + ^^n + l + 2^ #Vn} (130) 



where D = — \icx • V + j3m. From this we find the Cadzow equation of motion 



T^i = H D 7. 



2T 



(131) 
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where Hd = —ict ■ V + (3m and similarly for the conjugate field if)^- This is the 
required discretisation of the Dirac equation ( |127| ) . 

This equation has the same structure as that for the Grassmannian oscillator 
discussed in Paper I. To determine the regime for elliptic (particle like) behaviour, 
we first Fourier transform (|131|) , defining 

J d 3 xe-^ n (x) = # B (p) , (132) 

and then ( |131|) becomes 

(-3i + DT) ^ n+1 + (3^ + -DT) ^ = -4r/Z>„, (133) 

where D = a ■ p + /3m, which has the property 

DD = p • p + m 2 = £ 2 , (134) 

£7 being the definition of the particle's energy (keeping in mind there is no 'Hamilto- 
nian' in this theory). 

Now define the Dirac matrix 

3iDT + E 2 T 2 

v = , == 135 

^9E 2 T 2 + E 4 T 4 

and rescale the momentum fields via the rule 

= " n 4>n- (136) 

Then we find 

<f) n+1 + </> n _i = 2r/0„ (137) 

where 77 = -^===, which proves that oscillator behaviour is inherent in this system. 
Particle like solutions will occur for if < 1, from which we deduce the condition 
ET < v 7 ^. This upper limit is one half of the upper limit (|69|) established for the 
discretised Klein-Gordon equation. Exactly the same phenomenon occurs in the 
Grassmannian oscillator studied in Paper I. 

We shall leave the construction of conserved quantities such as the linear momen- 
tum and the Logan invariant built up from the ladder operators a (p) , b (p) , etc., as 
exercises. However, we will discuss here the charge and charge current densities. 

First, we rewrite the Cadzow equations of motion (|131| ) in the more useful form 

[-3t + TeT d ) ^ n+1 = -4Ti^V> n - (3i + TH* D ) ^ n _, 

(3z + H~ D T) = -^t%T - ^-1 (-3i + %T) , (138) 



where Hd = —ict ■ V + (3m, Hd = ict ■ V + /3m. 

Now the system function (|130|) is invariant to the infinitesimal global gauge trans- 
formation 

8^ n = i9j> n , 5^ + = -i9^+, (139) 
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so that there is a conserved charge given by 

Q n = J rf 3 x (| + y H2) VWi + ri+i (l - ^} • (140) 

It is easy to use the equations of motion (|138| ) to show that 

Q n = Q n -\ (141) 
The charge and current densities are then given by 

Pn = i>n (| + T Hd) ^n+l + V>n+1 (| ~ ) 

jn = \ fan «V»n-l + ^n a ^n + ^n-l a ^n} • ( M 2) 

These collapse to the standard densities in the limit T — > 0, and satisfy the discrete 
time equation of continuity 

Pn "Z"' 1 + v -j« = ° ( i43 ) 

as required. 

This completes our discussion of the free Dirac equation. The coupling of the 
Dirac field to the electromagnetic field in a gauge invariant way follows exactly the 
same procedure as with the charged Klein-Gordon field discussed in §8, with virtual 
paths following the rules given by ( |114j) . Second quantisation of the free Dirac field 
and QED will be discussed in the fourth in this series. 



10 Conclusions 

In this paper we have extended the methods outlined in the first paper of this series to 
classical fields, including scalar, vector, and spinor fields. The construction of equa- 
tions of motion and invariants of the motion is straightforward and gauge invariance 
also can be readily built into the system function. 

It is clear that the discretisation process does separate out timelike and spacelike 
terms in most equations, so that some of the elegance of a fully Lorentz covariant 
theory is lost. This is an inevitability in any discretisation process. However, this 
disadvantage should be viewed alongside the potential advantages, which may be 
considerable. 

First, our approach is exact, in that we do not countenance approximations within 
the theory itself, only in applications. Equations of motion are derived exactly from 
well-specified system functions, and invariants of the motion are precisely that, i.e., 
exact invariants of the motion and not approximate invariants of the motion. Next, 
there is a natural scale introduced into the theory, namely the fundamental time 
interval T. This influences the dynamics of the fields and leads to an upper limit for 
particle state energies, where the definition of particle state requires field amplitudes 
to remain bounded in time (corresponding to the condition rj 2 < 1 in the harmonic 
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oscillator). This holds some promise as a regulariser in quantum field theory, where 
in conventional calculations ad hoc cutoffs have to be introduced routinely. 

One of the surprises discrete time mechanics has sprung is how well it can work. 
Although it is certainly not equivalent to continuous time mechanics, in almost every 
aspect of the latter the former can make a comparable statement. We do not have 
a Hamiltonian in discrete time mechanics, but we do have a system function which 
can be used to generate equations of the motion, invariants of the motion, and define 
quantisation, as discussed in Paper /. In the case of free fields we can readily construct 
a Logan invariant which substitutes for the Hamiltonian in a particularly useful way, 
closed allied with the particle content of the theory. 

This raises the important question: if indeed all continuous time dynamical the- 
ories can be simulated by discrete time analogues, given a small enough T, so that 
no empirical test could distinguish between their predictions, what would force us to 
pick one approach rather than the other? One theory is elegant but badly defined 
mathematically in many places, whilst the other looks clumsier but is perhaps better 
defined mathematically. 

Ultimately it may come down to a question of personal choice, which will be 
influenced by some intrinsic features of the approach. For instance, the construction 
of Feynman path integrals as infinite products of integrals cannot be carried out 
in the usual way in discrete time field theory, simply because in that theory there 
is no meaning to the taking of the limit T — > 0, other than to show some sort of 
consistency with conventional theory. Therefore, one of the serious objections to the 
conventional path integral formalism, that of being ill-defined, is greatly moderated 
in discrete time theory. There is something like a path integral, but it involves a 
finite number of time steps between initial and final times, as discussed in Paper /. 
In Paper III of this series we shall discuss the second quantisation of field theories 
using this approach. 
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